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SUMMARY 

An analysis is presented for obtaining approximate solutions of the 
similarity equations for three-dimensional laminar-boundary- layer flows 
over a flat surface under main-flow streamlines that are translates and 
representable as infinite series expansions. For the particular case of 
streamline shapes described by a power of the distance along the surface 
from the leading edge, relatively simple expressions are obtained for 
flow deflection at the boundary surface, limiting streamline shape, and 
shear stress at the surface. 


INTRODUCTION 

In recent years, a great deal of attention has been focused on theo- 
retical investigations of three-dimensional incompressible boundary-layer 
flows. One class of investigations, having application to internal flow 
problems in turbomachines and flow over wings at high altitudes, has been 
concerned with finding exact solutions of the boundary-layer equations 
when the boundary layer develops over a flat surface (e.g., refs. 1 to 8). 
To date, all exact solutions have been based on similarity-type boundary- 
layer analyses. The essence of this technique involves the reduction of 
the partial differential equations for the boundary layer to a system of 
ordinary differential equations. The solutions for the actual boundary- 
layer flow are then obtained from the solutions of the ordinary differ- 
ential equations. It is generally necessary, however, to employ numerical 
methods and high-speed computing equipment to obtain accurate solutions 
to the ordinary differential equations because of their complex nature. 

As this process is time-consuming and often laborious, it is of interest 
to determine whether or not approximate solutions to such equations might 
be readily obtained which would encompass a wide variety of boundary- 
layer flows. 

Certain steps in this direction are taken in reference 9. In partic- 
ular, approximate solutions are obtained for the ordinary differential 
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equations arising from the similarity solution of the boundary-layer 
equations corresponding to the mainstream flows over a flat plate charac- 
terized by 


II 

tD 

(la) 

W = (const.) x 11 

(lb) 


(All symbols are defined in appendix A. The orientation of the coordi- 
nate axes is shown in fig. 1). The method of solution involves finding 
an approximate solution of the equations corresponding to the region of 
boundary-layer flow near the plate surface and then finding a separate 
approximate solution in the region near the main stream. The two indi- 
vidual solutions are then joined at a suitable point between these two 
regions (see also ref. 10, which advocates this technique) . Comparison 
of the approximate solutions of reference 9 with complete solutions ob- 
tained by relaxation techniques shows fairly good agreement when the value 
of n in equation (lb) is large (e.g., n = 10). This corresponds to the 
physical case of mainstream flows that flow for a distance without appre- 
ciable turning and then turn and accelerate rapidly, generating strong 
crossflows in the boundary layer. 

The present investigation also considers approximate solutions of the 
boundary-layer equations for the case of mainstream flows defined by equa- 
tions (l) . The reasons are twofold. First, although the approximate 
solution follows the procedure (advocated in refs. 9 and 10) of finding 
separate solutions of governing equations near the plate surface and near 
the main stream, the method of solution for the flow region near the plate 
surface differs from that given in reference 9 and hence provides a tech- 
nique of solution that is of interest in its own right. The method pre- 
sented herein is designed to give a more accurate approximation in the 
region near the plate surface and should yield improved estimates of shear 
stress and flow deflection at the surface. 

Secondly, reference 5 has shown that boundary-layer solutions for 
particular values of n in equation (lb) can be linearly combined to 
yield solutions of boundary-layer flows arising from mainstream flows 
given by 


o 

!3 

II 

13 

(2a) 

m 

w = 22 

(2b) 

n=0 


Because of computing-machine and time limitations, reference 5 obtains 
solutions only for a range of values of n from 1 to 10. This restricts 
the analysis to consideration of main-flow streamline shapes that can be 
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approximated lay eleventh-degree polynomials in x. Provided accurate 
approximate solutions can he obtained for values of n > 10 , the analysis 
of reference 5 could correspondingly he extended to include the calcula- 
tion of boundary layers for main-flow streamline shapes requiring approx- 
imating polynomials of degree higher than 11. For example, such approxi- 
mating polynomials would he required for circular-arc flows having more 
than 60° turning. 

The results of the present investigation and the approximate solu- 
tions of reference 9 are both compared with the exact solutions of refer- 
ence 5, and the nature of both approximate solutions is discussed. 

The work contained in this report is also included as part of a 
larger study presented in a doctoral thesis (ref. 11). Another portion 
of this thesis is presented in reference 12. 


ANALYSIS 

The solution for the boundary- layer velocity components for main 
flows over a flat surface defined by equations (2) is known to be (ref. 
5 ) 


u = UqF* (t|) 

( 3 ) 

1 (IF' - F) 

( 4 ) 

in 


= £ a^Pnfr) 

( 5 ) 


n=0 

where F(t]) is the well-known Blasius function satisfying the equation 


where 



= 0 



and P n satisfies the equation 


P" + 
r n T 


r n 


nF‘P n + n = 0 


( 6 ) 


( 7 ) 
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The boundary conditions on equations (6) and (7) are 

F(0) = F'(0) = P n (0) = 0 


(8) 


and 


lim F* (q) = lim P n (q) a 1 


( 9 ) 


T ) ’ > OD 


TH- 


Solutions to equation (7) (flows described by eqs. (2)) have been obtained 
for 0 < n < 10. Approximate solutions of equation (7) are now sought for 
all n > 10. To accomplish this, an approximate solution is first devel- 
oped for large values of q (i.e., in the vicinity of the main stream). 
Solutions are then found for small values of q (i.e., in the vicinity of 
the bounding surface). These two solutions are then matched at a suitable 
value of ri . 


An approximate solution of equation (7) for large t) and large val- 
ues of n is discussed in reference 13 and will be used in the present 
investigation. Basically, the solution is obtained as follows: For large 

values of T), it can be shown that P£ and P^ are necessarily small. 

Hence, if n is large, equation (7) can be written as F'P n * 1; and, 
therefore, the approximate solution for large tj and n can be given by 


Essentially the same result was obtained in reference 9 by considering 
the upper part of the boundary layer to be a "nonviscous" region and find- 
ing a solution of the boundary-layer equations with the viscosity term 
neglected. References 9 and 13 also point out that this approximate solu- 
tion applies over progressively wider ranges of values of q for progres- 
sively larger values of n. 


In order to obtain an approximate solution of equation (7) for small 
values of q, an adaptation of an analysis used in reference 14 to calcu- 
late skin-friction coefficients will be employed. In this regard, a new 
function «^(q) is defined, where ^ is related to P n by the equation 


Approximate Solution for Large q 



( 10 ) 


Approximate Solution for Small q 



( 11 ) 


Substituting equation (ll) into equation (7) and employing equation (6) 
give 
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.F”F* + ^2F" + - n(F') 2 Jk = -n (12) 

From the series expansion for the Blasius function F (see ref. 15, p. 
104) it is now assumed that F is adequately represented hy 

F = | (XT) 2 (a = F”(0) = 0.33206) 

in the neighborhood of T] = 0. Furthermore, examination of the second 
coefficient in equation (12) shows that FF'/2 is of the order of tj3 
relative to 2F"; and, hence, it is further assumed that in the neighbor- 
hood of t) = 0 the coefficient can be replaced by 2F” * 2a. Equation 
(12) can then be written as 



naqJk 


-n 

a,T) 


(13) 


At present, these approximations will be assumed sufficiently accurate 
for the purpose of the investigation. Later, results obtained from equa- 
tion (13) are compared with exact solutions in order to establish the 
validity of the assumptions. 


The solution of equation (13) can now be obtained analytically as 
follows . Let 

= Hn(0 (14a) 


where 


i 


2 i U 

= — n 5 a 5 T)2 


(14b) 


Substitution of equations (14) into equation (13) then yields 



(15) 


Consider the homogeneous complementary equation formed from equation (15) 


ff n,c + % £ - £ 2 + l) % c = 0 (16) 

Equation (16) is a form of Bessel's equation and has the general solution 

®n, c = ^n^l + ®n^l (17 ) 
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where Ii and Ki are Bessel functions of imaginary arguments (see ref. 

3 3 

16 for definitions and discussion). Tables of these functions are given 
in reference 17 . 


If the solution ^ c of the complementary equation (16) is known, 

the solution of equation (15) can he obtained by the method of variation 
of parameters. An analysis using this method is presented in appendix 
B, and the result is 


a _4 

= naa s 



Ii d£ 


Ii 

3 


f 

J o 


Ki d£ 

3 




V* 


(18) 


Hence, from equations (18), (14), and (ll), the solution for P n is 


PjH) = F* (h) t 



Ii d£ 


Ii 

3 



Ki d£ + A n Ii + BnKi 

(19 ; ) 


It is now recalled that, for regions near the wall, F(t]) was assumed 
to be adequately represented by the first term of its series expansion; 
that is, 

F 

2 


Using this approximation for F and using equation (14b), the expression 
for P n given by equation (19) becomes 



It is to be noted that equation (20) contains two arbitrary constants 
A n and B n . One of these constants can be determined by the boundary 

conditions P n (0) = 0. Setting t] equal to zero in equation (20) (equiv- 
alent to setting £ = 0) gives 
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where the series expansions for Ii(£) and Ki (£) (ref. 16) are used in 

3 3 

determining the equation. It follows at once from equation (2l) that 

Bn - 0. 


The determination of the remaining constant ^ depends upon the 

matching of the solutions for large and small r\ and is considered in 
the following section. 


Matching of Solutions 

The matching of the solutions for large and small t) consists in 
finding a value of rf where the values of P n and P^ for hoth cases 

can he brought into agreement. The equating of the functions and the 
derivatives gives a system of two simultaneous equations for the deter- 
mination of T] and the remaining arbitrary constant A n . 

The value of dP^dT] can be calculated from equation (20) by taking 
P n ^n d£ A i / 3\^ * , . 

liT = = -ar n3a3 (j) ^ 3 C 22 ) 


dt) 


Differentiating equation (20) with respect to £ gives 


dP E 

“IT 


(| r 1 K, + Kl)£ I, dt 


(i rl h + 1 i)fo K i d ^ A »(l rl V 1 ') 


Now, from reference 16 (p. 79), the following relations are known: 


(23) 


i t,' 1 Ki + Ki = -K_2 


(24) 


i r 1 ii + ii = i_2 

3 n n fl 


(25) 


Successive substitution of equations (24) and (25) into equation 
(23), and equation (23) into equation (22), gives 
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— = -a 3 ^ 


(i) 4 “ (■■• l 


Ki &£ + K_2 


f 


I* ^ - A n I 


-*) 


(26) 


The corresponding value of dP^/dri for the solution of the upper 

region of the boundary layer is determined from equation (10). Differ- 
entiating equation (10) with respect to q gives 


dR 


n 

dq 


(F’)‘ 


p" 


(27) 


At this point, it is assumed that the equations can be matched in a re- 
gion where F is reasonably approximated by F = ar\^f2. It will be 
shown later that this assumption will become increasingly accurate for 
larger values of n. Using this approximation for F and equating equa- 
tions "(10) and (20) and equations (26) and (27) give the following system 
of simultaneous equations for A n and the matching point q. (Actually, 

the equations are given as functions of £ and a matching value q 
can then be computed by means of the relation between £ and q given 
by eq. (14b).) 


(l) * ? ( K * 


H - 


r 


H 1 K i d E + ^ h 


) 


_1_ 

F* 


at] b 




2 

n3a 3 


(28) 


a U (|) ” C* (l-l jf d5 + K _2 l 


Ij. dt - A„ I 




F" 1 /s\~* A _i 

n8as 


_F" 

IF ) 2 


(29) 


From equation (28) it follows that 

/t 


r 1 + u 


An - 


r r 

I Ki - K, Ii 

Jo 3 3 Jo 3 


A ^ 

3 


II 

3 


(30) 
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Substituting equation (30) into equation (29) gives 

x -i | K i d ? + K -§£ H v - Tjjr 1 + I i| 5 k 4 a? - K i * 


1 j‘ § 


= 4 r 2 


(31) 


The solution of equation (3l) , which gives the value of £ at which 
the two solutions are matched, is outlined as follows: Multiplying equa- 
tion (3l) through by Ii and simplifying give 


% 


(v-s + V-s) ' l h « - E' 1 x -s - ! ?' 2 ij (32) 


Now, IiK_g + Kil_§ = (ref. 16, p. 80) . Therefore, equation (32) 
becomes 

% 


f 

Jo 


H d£ - 1-1 + f r 1 Ij 


(33) 


Numerical solution of equation (33) gives 

t, = X = 2.96 

It follows that q is given by 


(34) 


n = (£)' 


(l) a ' v 


- 3.902 

^ a. 

n 3 


(35) 


From the value of £ given in equation (34), it follows from equation 
(30) that the value of A^ is 


An = 1.811 


(36) 


Calculation of P n (ii) Curves 

From equations (20), (36), and the fact that % s 0, PnOl) can be 
expressed 
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P n (il) = Q(C) 


(37) 


where 


Q(E) - “ 1 (!)’ ( K i .f h ^ - h 


d£ + 1.811 I 


i) (“) 


Figure 2 shows a plot of Q(^) over the range 0<t;<£=2.96. In 
order to obtain a plot of P n (Tj) for a given value of n, it is necessary 

merely to select values of t], determine the corresponding values of £ 
from equation (14b), find Q(g) from figure 2, and then determine P n (rj) 
from equation (37 ) . 

A plot of the variation of tj with £ is presented in figure 3 for 
values of n equal to 10, 15, 20, and 50. Using these values of n and 
figures 2 and 3, curves for P n (Tj) were calculated and are presented in 
figure 4. 

As a check on the validity of the approximate solution, values of 
Pio(n) taten from reference 5 are plotted in figure 4. As the plot in- 
dicates, the values of the approximate solution fall slightly below the 
exact values for the "middle range" of T} values but agree very well 
for large T) and for tj near zero. 


Values of shear stress and limiting flow deflection near the bounding 
surface are of major interest in the analysis of three-dimensional 
boundary-layer flows. For the present case, both of these quantities can 
be shown to depend upon Pq( 0) (see ref. 5). In the following section, 

therefore, an expression for P^(0) is developed from equation (26), and 

values obtained are compared with exact values from reference 5 and the 
approximate solution of reference 9. 

Determination of P^(o). - From reference 16, the following expres- 
sions are noted; 


Flow Behavior Near Bounding Surface 



2 


+ higher -order terms 
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Hence, 


Employing 


Therefore, 


* (!i 


Ki = — + higher-order terms 

13 V3 r /2\ 

\3/ 


. (H 


K 2 = — 1 — 3 /, { + higher-order terms 

5 vsTgj 


M 

t 


1 1 _ — ^ ( + higher -order terms 

3 r /4\ 


I_2 

3 


K.2 

3 


fH 

o 

o 


Ki d£ = 


V3 « 

W 7 ® 


-f higher-order terms 


(i) a%/5 * * 

I* <1? - y L\ 7 1 V P + 


— Tj-r ^ 3 + higher -order terms 


4r (s) r (s) 


(39) and (40 ) in equation (26) gives 

n\-$ 


% - «**• (f 


"IU 


+ higher -order terms 


at the wall where £ = 0, 
fdP ' ~ -4 


/dP n \ a-int Aj, , 

= 2 - = 1.354 n 3 

V dTl A=o r (|)>^3 


(39) 


(40) 


( 41 ) 


The following table gives a comparison between the values calculated 
from equation (4l) and the exact values computed in reference 5: 
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n 

Pn(0)> exact 

values from 
ref. 5 

p A(o), 

eq. (41) 

i 

1.418 

1.354 

2 

2.197 

2.149 

3 

2.858 

2.816 

4 

3.450 

3.412 

5 

3.995 

3.959 

6 

4.506 

4.470 

7 

4.989 

4.954 

8 

5.449 

5.416 

9 

5.891 

5.858 

10 

6.317 

6.284 


It can be seen from the table that, for n - 8, 9, and 10, there is 
an error of approximately l/2 percent between the approximate solution 
i or P£(0) and the values presented in reference 5. 


, . A ^ point i-t is of interest to compare the values obtained for 

P n (0) in the present analysis with those obtained by the approximation 

technique in reference 9 . The method employed in the reference is based 
on expressing P n (ti) as a cubic polynomial in T) for small values of tj 

and as P n = 1 /f'(t]) for large values. The two solutions are matched in 

a manner similar to that presented in the present investigation. However 
as would be expected, the matching point tj is not the same for the two 
analyses. In reference 9 the matching point is given by 





It might be noted that, for very large values of n in equation (42), 
tion 3 (35)7* Th±S ±S t0 bS COmpaxed vith ^ " 3 - 9 02 n“* given in equa- 


p, (0 )f e “ aly8ls of reference 9 yields the following expression for 





( 43 ) 
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Figure 5 presents a plot of the values of P^(o) calculated from equation 
(41) and from equation (43). The ten exact values of P^(0) obtained in 

reference 5 are also plotted to show the agreement with the approximate 
solutions. The better agreement shown by the solution given in equation 
(41) may be accounted for as follows. First of all, the coefficient of 
A n in equation (41) is based on an exact solution of equation (13) . 

This is contrasted with the cubic polynomial representation in the neigh- 
borhood of q = 0 presented in reference 9. However, this must be 
weighed against the neglecting of the term FF"/2 in the second coeffi- 
cient of equation (12). For small q this approximation is reasonable. 
Secondly, A^ in equation (4l) depends upon the value of tj, or equiva- 
lently and it might, therefore, be inferred that the value of P^(0) 

would be critically influenced by the solution for large values of t] 
and the corresponding matching of the solutions for large and small tj. 
However, figure 6, which presents a plot of A n against £, shows that 

the value of A n is relatively insensitive to changes in for a wide 
range of values. Hence, it would appear that the determination of p^ (0 ) 

in the present analysis is principally affected by the accuracy of the 
technique used to evaluate the functions in the neighborhood of q = 0. 

Flow deflection at surface . - Assume that the mainstream velocity 
components are defined by equations (2). Then, if 7 designates the local 
angle of flow deflection at a point on the bounding surface of the flow, 
it follows from equations (3) and (5) that 


tan 7 = lim — 

T*0 


= lim 
T]->0 


£6 F'(tD 


v 4*^(0) , . 

~ ifo F” (0) ( 44 ) 

where and where L'Hospital's rule is used to evaluate the 

indeterminate form P n (0)/F‘(0). 


As equation (4l) has been shown to yield a reasonable approximation 
to Pn(0) for all n, equation (44) can be written as 


tan 7 £ 



(45) 
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where 



1.354 


* 

an «- 

F”"(0 ) = 4,08 


For the special case 


W » a n x 11 

u - u 0 


equation (45) becomes 


tan y = b^x^n^ 


(46) 


(47) 


Limiting streamline . - The equation for the "limiting streamline” 
(line of limiting flow deflection, ref. 5) at the surface can be deter- 
mined directly from the relation 


dz 


dx 


■ £ “*■-* 


where z^ is the z-coordinate of the limiting streamline. 
Solution of equation (48) gives 


(48) 


z , * V b * n+i . n 1 

^ > b n x — i 


+ constant 


(49) 


as the approximating equation of the limiting streamline. 

For flows defined by equation (46), equation (49) becomes 


* ntl (dr) 


+ constant 


Now, the main-flow streamlines for this case are given by 


a£x n+1 

z = — =— + constant 

n + 1 


If the constant of integration is chosen equal to zero, the ratio of z^ 
to z is given by 
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~i - = -y n$ = 4.08 n$ 
Z a n 


( 50 ) 


Hence, the ordinates of the boundary-layer limiting streamline can be 
determined at once from the ordinates of the main-flow streamline through 

the origin merely by multiplying the latter by the scale factor 4.08 n$. 
The entire system of limiting and main-flow streamlines can then be ob- 
tained by translation parallel to the leading edge. 


Shear stress at surface . - The shear stress at the bounding surface 
is given by 


T 0 “ |/( T 0,x) 2 + ( T 0,z) 2 

where (using eqs. ( 3 ) and (5) and the definition of rj) 


0,x = (“ s) y=0 - W'(°) Vi 


(51) 


(52) 


/ 

c °> z = ( M 3?) y=0 “ 11 § 4* n PA(o) V^ 


(53) 


It follows from equations (51), (52), (53), and (41) that t q can be 
expressed as 


t 0 « — = 1 / f 


/vx 


V F "( 0) + [n?0 a S xI1 d.354 n$)J 


For flows defined by equation (46), equation (54) becomes 


H.U n 2 t / 

= ~= ^F"(0)n-$ + (4) 2 (1.354 x*)2 


Vvx 

Therefore, for sufficiently large values of n, 

pU n ^ 1 

T 0 = -y/vx a n(-'- , 354)x I1 n 3 


(54) 
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CONCLUDING REMARKS 

A method has been developed for obtaining approximate solutions of 
the similarity equations for three-dimensional laminar -boundary-layer 
flows over a flat surface under main -flow streamlines that are translates. 
The approximate solution makes possible the analysis of flows having 
streamline shapes representable by polynomials of any degree. 

For the particular case of streamline shape described by z = ax 11 , 
relatively simple expressions can be found for flow deflection at the 
surface, limiting streamline shape, and shear stress at the surface. 

Comparison of the approximate solution with exact solutions obtained 
on high-speed computing equipment shows good agreement for values of 
n = 10. In particular, the approximate solution in the vicinity of the 
plate surface shows relatively good agreement with the exact solution for 
all values of n. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 23, 1956 
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APPENDIX A 
SYMBOLS 

A n 

constants of integration 


constants 

B n 

constants of integration 

*n 

constants 

F,F(t]) 

Blasius function (eq. (6)) 


function of similarity parameter (eq. (ll)) 


function of £ (eq. (14a)) 

I V ,K V 

Bessel functions of imaginary argument (eq. (17)) 

m,n 

constants 


function of similarity parameter (eq. (7)) 

Q(C) 

function of £ (eq. (38)) 

u,w 

components of mainstream velocity in x- and z-directions, 
respectively 

u o 

inlet velocity of main stream 

u,v,w 

components of "boundary- layer velocity in x,y, z directions , 
respectively 

x,y,z 

rectangular coordinates 

a 

constant, - F" (0) = 0.33206 

r 

gamma function 

r 

flow deflection angle at surface 

C 

function of similarity parameter t) (eq. (141)) 
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T| similarity parameter,, T) = y \/Uq/vx 

p coefficient of absolute viscosity 

v coefficient of kinematic viscosity 

Tq shear stress at wall 

Subscripts: 

c value of complementary function (eq. (16)) 

l limiting 

m constant 

n index number 

x,z x- and z-components, respectively 

v order of Bessel functions 

Superscripts: 

' differentiation 

* division by U (ref. 5) 

— value at matching point 
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APPENDIX B 

DERIVATION OF EQUATION (l8) 

Let Sj_ and s 2 be solutions of the homogeneous equation obtained 

from 


Y" + p(x)Y ' + q(x)Y = r(X) 


(Bl) 


by setting r equal to zero. Then, reference 18 shows that the general 
integral of (Bl) is 

C rs 2 rrs x 

Y = -si/ ~ dX + s 2 J clX 

where S = s^Sg - s^s^- 


Applying this technique to equation (15) and using I a and Ka as 
solutions of (16) give the following solution for H n : 



d£ + Ki 

3 3 



-1 

TIk 

5 


1 

3 



(B2) 


Reference 16 shows that 

IiKi - IiKi = -r 1 (B3) 

3 3 3 3 

Substituting (B3) into (B2) then gives 

= n$(§) a-^Ii/K, d£ + K./li dtj (B4) 

By expressing the integrals in equations (B4) as definite integrals with 
additive constants, equation (18) is obtained. 
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Figure 1. - Coordinate axes for flow over surface. 
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Figure 3 


Variation of similarity parameter 7? with variable £ 
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Figure 4. - Approximate solutions of P n ( t}) . 
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Figure 5. - Comparison of approximate and exact solutions for P'(0). 
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Figure 6. - Variation of A n with function £ at matching point. 
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